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Section 


ABSTRACT 


This  report  presents  a  novel  Kahnan  filtering  approach  to  the  suppression  of  narrowband 
interference  from  direct  sequence  spread  spectrum  communications  systems.  The  algo¬ 
rithm  is  based  on  the  digital  phase-locked  loop  Kalman  filter.  Because  the  interference  is 
assumed  to  be  much  stronger  than  either  the  signal  or  noise,  the  Kalman  filter  locks  onto 
a  function  related  to  the  interference.  The  net  result  is  an  estimate  of  the  phase  of  the 
interference  and  its  amplitude.  The  algorithm  is  characterized  through  computer  simu¬ 
lation  for  the  case  of  narrowband  Gaussian  noise  interference.  Examples  of  the  phase- 
and  envelope-tracking  capabilties  of  the  algorithm  are  presented,  followed  by  bit  error 
rate  curves  for  interference  bandwidths  ranging  from  0.5%  to  5.0%  of  the  chip  rate  of  the 
spread  spectrum  signal  .  The  report  concludes  with  three  adaptive  architectures.  The 
first  is  suitable  for  constant  envelope  interference;  the  second  is  a  more  general  structure; 
and  the  third  incorporates  a  decision-feedback  structure  accompanied  with  a  traiining 
sequence. 


RiisuM^: 

Ce  rapport  presente  une  nouvelle  fa^n  d’effectuer  im  filtrage  de  Kalman  pour  exciser 
de  I’interference  a  bande  etroite  d’un  systeme  de  communication  utilisant  des  spectres 
etciles  par  sequence  directe.  L’algorithme  est  base  sur  un  filtre  de  Kalman  a  boucle 
numerique  de  verrouillage  de  phase.  L’interference  est  sens^  etre  beaucoup  plus  faible 
que  le  signal  ou  le  bruit  de  fa^n  a  ce  que  le  filtre  de  Kalman  se  verrouille  sur  une 
function  liee  a  I’interference .  II  en  resulte  une  evaluation  de  la  phase  et  de  I’amplitude  de 
I’interference.  L’  algorithme  est  caracterise  grace  a  une  simulation  pax  ordinateur  pour  le 
cas  ou  I’interference  prend  la  forme  de  bruit  gaussien  a  bande  etroite.  On  pr^nte  d’abord 
des  exemples  de  poursuite  de  phase  et  d’enveloppe,  puis  des  courbes  du  taux  d’erreur 
des  bits  pour  des  bandes  d’interference  ayant  une  largeur  entre  0.5%  et  5.0%  du  debit 
numerique  du  signal  a  spectre  etale.  Le  rapport  se  termine  par  la  presentation  de  trois 
architectures  adaptives.  La  premiere  convient  aux  interferences  a  enveloppe  constante, 
la  deuxieme  est  une  structure  plus  generale  et  la  troisieme  incorpore  une  structure  de 
retroaction  decisionnelle  et  une  sequence  d’essai. 
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EXECUTIVE  SUMMARY 


This  report  presents  a  Kalman  filtering  algorithm  that  is  used  for  filtering  narrowband 
interferers  out  of  direct  sequence  spread  spectrum  signals.  These  signals  are  used  exten¬ 
sively  in  military  communication  systems.  The  technique  described  herein  applies  equally 
to  both  Electronic  Support  Measures  (ESM)  systems  and  direct  sequence  spread  spectrum 
communication  sytems.  In  the  former  application,  the  ESM  system  may  be  attempting  to 
intercept  the  spread  spectrum  signal,  but  the  narrowband  interference  may  be  hampering 
this  effort.  In  the  latter  application,  the  spread  spectrum  communication  system  may 
require  additional  assistance  to  suppress  the  interference.  Since  the  open  literature  has 
been  devoted  to  this  latter  case,  the  material  presented  here  focuses  on  this  application. 

One  of  the  attributes  of  direct  sequence  spread  spectrum  communication  systems 
is  their  ability  to  combat  interference  or  intentional  jamming  by  virtue  of  the  system’s 
processing  gain  inherent  in  the  spreading  and  despreading  process.  The  interference  can  be 
attenuated  by  a  factor  up  to  this  processing  gain.  In  some  cases  the  gain  is  insufficient  to 
effectively  suppress  the  interferer,  leading  to  a  significant  degradation  in  communications 
manifested  by  a  sudden  increase  in  bit  error  rate.  If  the  ratio  of  interference  bandwidth 
to  spread  spectnim  bandwidth  is  small,  the  interference  can  be  filtered  out  to  enhance 
system  performance.  However,  this  is  at  the  expense  of  introducing  some  distortion  onto 
the  signal.  This  process  of  filtering  is  sometimes  referred  to  as  interference  excision. 

The  Kalman  filtering  algorithm  is  based  on  the  digital  phase-locked  loop  Kalman 
filter.  Because  the  interference  is  assumed  to  be  much  stronger  than  either  the  signal  or 
noise,  the  Kalman  filter  locks  onto  the  interference.  The  net  result  is  an  estimate  of  the 
phase  and  amplitude  of  the  interference. 

The  algorithm  is  characterized  through  computer  simulation  for  the  case  of  nar¬ 
rowband  Gaussian  noise  interference.  Examples  of  the  phase-  and  envelope-tracking  ca- 
pabilties  of  the  algorithm  are  presented,  followed  by  bit  error  rate  curves  for  interference 
bandwidths  ranging  from  0.5%  to  5.0%  of  the  chip  rate  of  the  spread  spectrum  signal. 
The  report  concludes  with  three  adaptive  architectures.  The  first  is  suitable  for  constant 
envelope  interference;  the  second  is  a  more  general  structure;  and  the  third  incorporates  a 
decision-feedback  structure  used  in  conjunction  with  a  truning  sequence  to  stabilize  the 
filter. 

Military  applications  of  this  work  cannot  be  discussed  in  this  unclassified  re¬ 
port. 
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1.0  INTRODUCTION 


Direct  sequence  spread  spectrum  commimication  systems  have  an  inherent  processing 
gain  which  can  reduce  the  effects  of  jammers  or  unintentional  interference.  When  these 
intruding  signals  have  a  power  advantage  over  the  spread  spectium  system,  a  severe  degra¬ 
dation  in  communications  results.  However,  communications  can  be  enhanced  somewhat 
by  filtering  the  interference,  particularly  if  its  bandwidth  is  significantly  less  than  the 
bandwidth  of  the  spread  spectrum  signal. 

A  host  of  digital  filtering  algorithms  based  on  time-modelling  concepts  have 
been  developed  over  the  years  to  address  this  problem  [1,  2,  3].  The  received  signal, 
noise  and  interference  are  applied,  for  example,  to  a  filter  matched  to  a  chip  and  the 
output  is  sampled  at  the  chip  rate.  If  one  assumes  that  the  resulting  sampled  data  can 
be  modelled  as  an  autoregressive  (AR)  process,  then  the  AR  coefficients  or  their  related 
lattice  reflection  coefficients  can  be  determined  from,  say,  the  maximum  entropy  algorithm 
[4],  stochastic  gradient  algorithms  [5],  and  block  and  recursive  least  squares  algorithms 
[6,  7].  Estimates  of  these  coefficients  by  the  above  algorithms  lead  directly  to  transversal 
or  lattice  filter  structures,  both  of  which  act  as  whitening  filters  for  the  assumed  AR 
process.  This  approach  to  suppress  interference  (or  excision  as  it  is  sometimes  referred 
to)  is  possible  because  of  the  non-coherency  of  the  signal  and  noise  samples,  and  the  more 
coherent  interference  samples  resulting  from  the  latter’s  narrowband  nature. 

This  paper  presents  a  different  time-modelling  approach  to  interference  excision. 
A  function  related  to  the  interference  is  assumed  to  have  been  generated  by  a  second 
order  state-space  system  to  which  the  extended  Kalman  filter  equations  are  applied.  This 
results  in  what  has  been  termed  the  digital  phase-locked  loop  (DPLL)  [8,  9].  There  are 
certain  benefits  to  using  the  DPLL  over  its  analog  counterpart  [10].  First,  it  can  be  made 
adaptive,  responding  to  changing  interference  conditions;  and,  second,  being  digital,  other 
DPLL’s  can  be  switched  in  as  required  to  handle  more  than  one  interferer. 

The  outline  of  the  paper  is  as  follows.  Section  2.0  describes  the  spread  spec¬ 
trum  communication  model.  Section  3.1  describes  the  state-space  model  which  produces 
a  function  related  to  the  interference.  Section  3.2  presents  the  Kalman  algorithm  which 
is  used  to  estimate  the  phase  of  the  interference.  Another  model  [11]  is  also  presented, 
which  is  slightly  different  from  that  in  [8].  Section  4.0  describes  the  interference  estima¬ 
tor.  Section  5.0  presents  simulation  results  for  the  case  of  narrowband  Gaussian  noise 
corrupting  the  spread  spectrum  signal.  The  results  are  presented  in  terms  of  bit  error 
rate  for  various  parameter  values  of  the  excisor.  Finally,  the  results  suggest  a  strategy  to 
render  the  algorithm  adaptive,  which  is  discussed  in  Section  6.0. 


2.0  MODEL  OF  THE  PN  SPREAD-SPECTRUM  COMMUNICATION 

SYSTEM 


A  BPSK  PN  spread  spectnim  system  is  considered  here,  the  basic  elements  of  which  are 
shown  in  Fig.  1. 


pn 

Correlator 


Figure  1:  Spread  spectnim  communications  model. 

Certain  assumptions  are  made  in  what  follows.  These  include  distortionless 
transmission,  and  phase  and  code  synchronization. 

The  received  waveform  r(t),  consisting  of  a  spread  spectrum  signal,  additive 
white  Gaussian  noise,  and  narrowband  interference  is  applied  to  a  bandpass  filter  with 
the  transfer  function  whose  output  is  defined  as 

z(t)  =  s(t)  +  n(t)  +  *(t).  (1) 

The  bandpass  filter  for  the  application  considered  in  this  paper,  is  zissumed  to  be 

a  filter  matched  to  a  chip  and  centered  at  the  carrier  frequency  u>o  of  the  spread  spectrum 
signal. 

The  spread  spectrum  signal,  which  is  of  bandwidth  is  defined  as 

s(t)  =  a(t)cos(wot),  (2) 

where 

L 

=  ES CnjdnPit  -  jTc  -  nTi).  (3) 

«  jatl 

The  term  Cnj  in  Eq.  (3)  is  the  code  sequence  for  the  n*^  information  bit  d„,  Th  is  the 
reciprocal  of  the  bit  rate  Rt,,  p{t)  is  the  basic  chip  pulse  of  energy  Ec,  Te  is  the  reciprocal 
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of  the  chip  rate  Rc,  and  L  is  the  number  of  chips  per  bit,  also  known  cs  the  processing 
gain. 

The  noise  n(<)  is  Gaussian  and  has  a  power  spectral  density 

5.(/)  =  w 

where  iVo/2  is  the  power  spectral  density  of  the  assumed  white  Gaussian  noise  from  the 
channel.  The  interference  is  generally  defined  as 

t(f)  =  /(f)cos[wbf  +  ^(i)l  (5) 


and  is  of  bandwidth  Bi  «  B,,. 

Referring  to  Fig.  1,  the  output  z{t)  of  the  bandpass  filter  /f6p(/)  is  bandpass 
sampled  and  applied  to  a  bandpass  limiter  and  interference  estimator. 

Consider  the  bandpass  sampler  first.  The  analog  signal  z{t)  from  Eq.  (1)  is 
sampled  at  f,  =  2Rc  {mf,  =  wo/25r  +  Rc  for  some  integer  m),  where  Rc  is  the  chip  rate. 
The  resultant  sampled  signal  is,  therefore, 

z(ib)  =  s(ifc)  +  n(/:) +  *(*),  (6) 

where  s(ifc)  consists  of  the  sequence  {. . . ,  0,  (-1)*-Ec>  0,  (— •  •  •}  (recall  that  coher¬ 
ent  bandpass  sampling  has  been  assumed),  n{k)  are  uncorrelated  Gaussian  noise  samples 
of  variance  <7^  =  EciNof2),  and  i{k)  is  the  sampled  version  of  Eq.  (5).  The  samples  z{k) 
are  applied  to  the  interference  estimator  and  interference  canceller. 

Consider  now  the  branch  containing  the  limiter.  Here,  z{k)  is  applied  to  a 
bandpass  limiter.  The  input  to  the  limiter  is  redefined  as 

z(k)  =  yj[lik)  +  nl(fc)  +  cosfu;*,*  +  +  'A.(*)]  (V 


where 


^^(Jb)  =  arctan 


/  n'M-haW  \ 

[l(k)  +  n'i(k)  +  a'i(k)J 


(8) 


is  a  noiselike  phase  fluctuation  on  the  interferer’s  phase  0(k),  and  is  due  to  the  noise  and 
spread  spectrum  signal.  The  terms  n[(k),  n^lk),  and  a3(^)  are  in-phase  and  quadra¬ 
ture  components  of  the  noise  and  spread  spectrum  signal  with  respect  to  the  interference 
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phase  6{k).  The  output  of  the  bandpass  limiter  is  [12] 

bp{k)  =  ^  cosfwoib  +  9(k)  +  ^,ik)],  (9) 

where  A'  is  the  limiter’s  output  level.  This  signal  is  redefined  as  (letting  A'  =  \/2ir/4  for 
convenience) 

bp{k)  =  V^cosfo^A;  +  d(k)  +  ^x(k)].  (10) 

It  should  be  noted  that  for  large  interference-to-noise  ratios  in  which  the  interference  is  of 
constant  envelope,  <l>x{k)  in  Eq.  (9)  is  approximately  Gaussian  [12].  The  sampled  signal 
in  Eq.  (10)  is  what  is  processed  by  the  Kalman  filter,  which  estimates  the  phase  6{k)  of 
the  interference. 

3.0  KALMAN  FILTER  ALGORITHM 

3.1  STATE  SPACE  MODEL  OF  THE  INTERFERENCE 

The  development  here  follows  that  in  [8].  The  analog  interference  model  will  be  discussed 
first,  followed  by  its  digital  counterpart. 

A  block  diagram  of  the  analog  state  space  model  which  generates  the  analog 
version  of  a  function  related  to  6p(<)  in  Eq.  (9)  is  shown  in  Fig.  2.  The  state  space 
equations  are  defined  as 

xi(t)  1  _  r  0  dir  xi(<)  ]  [  0 

h(t)  J  [  0  -a  ]  [  n(()  J  [  vTC  J  *  *  * 

where  Xi{t)  and  X2(t)  are  the  phase  function  and  modulating  signal  state  variables  of  the 
interference,  respectively.  In  compact  form  Eq.  (11)  becomes 

i(t)  =  Ax{t)  +  gw{t).  (12) 

The  interferer’s  modulating  signal  X2(t)  is  assumed  to  be  ge  *  Aed  by  appl}ring  Gaussian 
white  noise,  ti;(t),  of  unit  variance  to  an  amplifier  with  gun  y/KZ  followed  by  a  first  order 
low  pass  filter  of  bandwidth  a  =  2xa/.  The  steady  state  variance  of  the  coloured  noise 
at  the  output  of  the  lowpass  filter  is  unity  if  K^,  =  2a.  The  output  is  multiplied  by  the 
frequency  deviation  constant,  d,  and  then  integrated  to  yield  the  phase  function 

Xi{t)  =  d  f  X2(t)  dr.  (13) 

Jo 
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An  FM  type  of  signal  (interferer)  is  generated  by  phase  modulating  a  carrier,  ujq,  which 
has  been  set  equal  to  the  carrier  of  the  spread  spectrum  signal,  since  any  offsets  from  u>o 
can  be  incorporated  into  yielding 

y(t)  =  V5co8[wft< +  ari(t)].  (14) 

This  signal  is  corrupted  by  white  Gaussian  noise  v{t)  with  power  spectral  density  N/2 
generally  different  from  No/2  defined  earlier  for  the  communications  model  in  Fig.  1.  The 
result  is  an  FM  signal  related  to  the  interference  in  Eq.  (9)  except  for  the  faM:t  that  the 
noise  term  is  additive  here,  i.e,, 

bpity  =  yit)  +  vit).  (15) 

For  large  interference-to-signal-plus-noise  conditions,  the  final  effect  of  the  additive  noise 
in  estimating  the  phase  xi(t)  in  Eq.  (14)  will  be  similar  to  the  effect  of  in  Eq.  (9). 

Given  the  observation  process  bp\t)  from  Eq.  (15),  the  objective  is  to  estimate 
the  state  Xi(t)  and  the  related  state  X2(t),  as  well  as  y{t)  in  Eq.  (14).  One  approximate 
technique  reserved  for  nonlinear  estimation  problems  [13]  is  the  extended  Kalman  filter 
algorithm. 


v(t) 


Figure  2:  Block  diagram  of  the  state  space  model. 


Before  the  2tlgorithm  is  presented,  the  discrete  forms  of  Eqs.  (11)  and  (12)  will 
be  stated,  namely, 

»(Ar  +  1)  =  <^aj(A:)  +  «;(*),  (16) 

where  k  represents  discrete  time, 

1  [l  <f>u' 

0  0  ^2 
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(17) 


and 


(18) 


w{k)  =  /  <f>gw{u)  du. 
JkT, 


In  Ek}.  (17),  ^  =  d/a  and  is  termed  the  bandwidth  expansion  ratio  in  units  of  volts~^ 
[14].  In  Eq.  (18),  w{k)  is  a  stationary  zero-mean  white  Gaussian  vector  sequence  whose 
covariance,  (i.e.,  ^{m(A:)tD(A:)‘}),  is  given  by 


f(k+iYr, 

Vw=  (f>gg*4>du.  (19) 

JkT, 

Finally,  Eq.  (14)  becomes, 

h'{k)  =  y(fc)  +  v{k),  (20) 

where  y{k)  is  the  sampled  version  of  Eq.  (14).  It  has  been  assumed  that  bj/{t)  from 
Eq.  (15)  has  been  applied  to  a  bandpass  filter  of  bandwidth  B,  whose  output  has  been 
sampled  at  a  rate  so  as  to  yield  uncorrelated  noise  samples  v{k)  of  variance  =  B{N/2) 
[9]. 

The  noise  covariance  matrix,  V^,,  is  obtcuned  by  substituting  Eq.  (17)  and  g 
from  Ekj.  (12)  into  Eq.  (19)  and  integrating  each  term,  yielding 


^  (47r/7  -  3  + 

2a  ^  (l  -  1  - 


(21) 


where  7  =  2ir/{aT,),  which  is  the  ratio  of  the  sampling  rate  f,  =  1/T,  to  the  interferer’s 
modulating  signal  bandwidth  aj.  The  discrete  form  of  the  interference  model  is  shown  in 
Fig.  3. 


Figure  3:  Discrete  form  of  the  interference  model. 
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3.2  EXTENDED  KALMAN  FILTER  EQUATIONS 


The  extended  Kalman  filter  ^nations  [13]  can  now  be  applied  to  the  discrete  model.  The 
resulting  algorithm  is  as  follows: 


Initial  Conditions: 


Do  for  fc  =  0  to  T: 


e-Hib|ib-l)  = 


x(-l)  =  E{xi-1)}  = 
Vi(-l)  =  E{xi-l)x\-l)} 


e{k\k-l)=[  1  0  ]0i(ib-l) 
Vi{k\k  -  1)  =  <l>Vi{k  - 
_ 1  —  co8[2u^<;  +  20{k\k  —  1)] 


<^2+[l  0  ]  -  1) 


1  0 
0  0 


Vi(ib|ib-l)e-nifc|Ar-l) 


K{k)  =  -V,{k) 


v^sin[u^fc  +  §{k\k  —  1)] 


x{k)  =  0x(fc  —  1)  +  lir(fc)[5p'(fc)  —  >/2cos[u^fc  +  9{k\k  —  1)]. 


(22) 

(23) 

(24) 

(25) 

(26) 


(27) 

(28) 


In  the  algorithm  x{k)  is  defined  as  the  estimate  of  the  true  state  vector  x{k),  x{k) 
is  defined  as  the  post  error  x{k)  -  x{k),  K{k)  is  the  Kalman  gain  vector,  Vi^k)  = 
E{x{k)x*{k)]  is  the  post  error  covariance  and  Vi{k\k  -  1)  =  E{x{k\k  -  l)i‘(ib|Jfc  -  1)}. 
Furthermore,  Elq.  (26)  is  known  as  the  “coupled  equation”,  since  the  estimate  of  the 
interferer’s  phase  is  contained  in  it.  However,  this  equation  can  be  further  approximated 
as  noted  in  [8].  The  approximation  is 

It  has  been  shown  in  [8]  that  the  K^Jman  algorithm  results  in  a  digital  PLL.  The 
state  space  diagram  of  the  DPLL  is  illustrated  in  Fig.  4.  Observe  that  the  phase  estimate 
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6{k\k  —  1)  is  applied  to  a  phase  modulator  of  frequency  resulting  in  a  “VCO”  signal 
defined  as 

VCO(k)  =  - v^sin(a;oib  +  e{k\k  -  1)).  (30) 

Shifting  this  signal  by  — x/2  radians  results  in  a  filtered  estimate  of  y(ib)  in  Eq.  (20). 
Finally,  by  rearranging  the  system  in  Fig.  4  into  a  transfer  function  model  at  steady 
state,  one  obtsuns  the  system  illustrated  in  Fig.  5  where 


^  2^„V*..(oo)  +  Vj„(oo) 

(31) 

.  _  ^12l4„(oo)  +  (00)^2 

2^12V*„(oo)  +  V*„(oo) 

(32) 

a  =  (^22- 

(33) 

The  DPLL  consists  of  a  gain  C,  followed  by  the  loop  filter,  and  integrator.  These  equations 
indicate  that  the  gain  and  loop  filter  pau'ameters  (pole  and  zero)  are  linked  to  the  steady 
state  phase  variance  14,j(oo),  the  phase/modulating  signal  covariance  Vi,j(oo),  the  state 
transition  parameters  and  ^22,  and  the  variance  of  the  observation  noise  <t*.  Steady 
state  performance  ultimately  depends  on  the  choices  of  d,  and  T,  of  the  model. 


Figure  4:  State  spiu:e  representation  of  the  DPLL. 


The  loop  filter  shown  in  Fig.  5  will  lead  to  a  steady  state  error  if  the  phase  xi(jb) 
is  a  ramp.  This  loop  filter  can  be  changed  to  one  consisting  of  proportional  plus  integral 
control  by  modifying  the  interference  model  [11].  The  modification  involves  replacing  the 
low  pass  filter  by  an  integrator  ( i.e.,  by  setting  a  =  0  in  A  in  Eq.  (12)).  With  this  change. 
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Figure  5:  Transfer  function  model  representation  of  the  DPLL. 


the  state  transition  matrix  and  covariance  matrix  become,  respectively, 


1  dr/ 

1  <f>12 

0  1 

0  <f>22 

(34) 


and 


V^  =  K^ 


<PT^  dT^ 

dr*  T, 


(35) 


Using  these  terms  in  the  Kalman  equations  results  in  a  slightly  modified  algorithm.  This 
is  the  algorithm  that  has  been  used  in  the  simulations  in  Section  5.0.  The  transfer  function 
parameters  in  Eqs.  (31)  and  (32  )  remain  imchanged,  but  a  =  1  in  Eq.  (33). 


4.0  THE  INTERFERENCE  ESTIMATOR 

The  Kalman  algorithm  produces  a  filtered  estimate  of  y{k)  in  Eq.  (20),  which  can  be 
used  as  a  basis  for  estimating  a  sampled  version  of  the  the  amplitude  of  t(f)  (i.e.,  I{k)  in 
Eq.  (5)).  The  interference  estimator  is  illustrated  in  Fig.  6.  The  output  of  the  multiplier 
is  a  baseband  term,  and  is  (using  Eq.  (6)), 

/'(t)  =  -^{/(i)ca.(tf(t)-*(i|t-l)) 

+  n{k)  co6[ijfok  +  —  1)]  +  <*(^)  cos[(?(jb|Jb  —  1)] 

+  o(A:)  co8[2woA:  +  ^(fc|t  -  1)]}.  (36) 

Equation  (36)  consists  of  four  terms:  the  first  term  is  related  to  the  desired 
amplitude  of  the  interference;  the  second  is  approximately  Gaussian  baseband  noise  [14]; 
and  the  third  and  fourth  terms  are  noise-like  terms  emanating  from  the  spread  spectrum 
signal.  The  fourth  term  has  been  ignored  in  the  simulations,  since  it  has  been  assumed 
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Figure  6:  Block  diagram  of  the  interference  estimator. 

that  it  contributes  insignificant  aliased  noise  into  the  baseband  region  as  a  result  of  the 
sampling  rate  conditions  discussed  in  Section  2.0.  The  term  r{k)  is  filtered  by  the  low  pass 
filter  of  bandwidth  Blpfi  resulting  in  the  estimate  of  the  interference  envelope,  i{k)l\/2. 
Combining  this  with  y{k)  shown  in  Fig.  6  yields  the  estimate  of  the  interference,  i.e., 

i{k)  =  i{k)  cosfcJbAr  +  9{k\k  —  1)],  (37) 

which  is  subtracted  from  z{k)  to  produce  the  error  signal  e(k)  as  illustrated  in  Fig.  1. 
Because  the  input  was  coherently  bandpass  sampled  at  twice  the  chip  rate,  this  error 
signal  tmdergoes  decimation  and  alternate  sign  reversal  first,  followed  by  correlation  with 
the  local  PN  code  and  bit  detection.  Several  examples  of  the  performance  of  the  estimator 
will  be  presented  next  for  the  case  of  narrowband  Gaussian  noise  interference. 


5.0  SIMULATION  RESULTS 


This  section  describes  the  simulation  conditions  of  tests  which  were  conducted  and  presents 
examples  of  the  performance  of  the  exdsor,  including  phase  tracking,  envelope  estimation, 
and  bit  error  rates.  The  parameters  used  in  the  simulations  are  listed  in  Table  1. 

The  modified  Kalman  algorithm,  as  alluded  to  by  Eqs.  (34)  and  (35),  has  been 
used  in  the  simulations.  The  initial  conditions  were  as  follows: 


*(-l)  =  E{xi-1)}  = 


(38) 


0.1  0.0 

0.0  0.1 


(39) 


The  choice  for  the  initial  values  in  Eq.  (39)  is  as  follows.  If  the  input  to  an  integrator  is 
zero  mean  white  Gaussian  noise  of  some  variance,  the  output  variance  will  increase  linearly 
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Table  1:  Simulation  parameters.. 


Signal 

Interference 

Kalman  Filter 

i  =  20 

E^INq  =  10  dB 
Ec  =  \  unit 
i2e  =  1  Hz 

Bi  =  0.5%  to  5%  of  Rc 

W,*  =  (jjQ 

I/S  t(per  chip)  =  20  dB 

a  =  2x(.0003125)  rad./sec. 

=  2a 

d{k)  =  t 

N/2  =  257,  102,  and  32.4  units 

*  Carrier  frequency  of  interference, 
t  Interference-to-Signal  ratio. 

^  A  time-varying  frequency  deviation  constant  is  used  in  one  example. 


with  time,  starting  from  an  initial  value  of  zero.  If  the  output  of  the  first  integrator  is 
applied  to  a  second  integrator,  then  the  variance  of  the  output  of  the  second  integrator  will 
increase  quadratically  with  time,  also  starting  with  cm  intial  value  of  zero.  To  start  the 
simulation  with  a  positive  definite  covariance  matrix,  small  positive  values  (i.e.,  0.1)  were 
inserted  on  the  diagnal  of  Eq.  (39).  It  should  be  noted  too  that  the  impact  of  the  initial 
starting  value  will  be  transitory,  eventually  settling  down  to  a  steady  state  value.  Finally, 
the  bandwidths  for  the  low  pass  filter  Blpf  shown  in  Fig.  6  and  used  in  determining  the 
envelope  of  the  interference,  ranged  from  0.1  Hz  to  0.4  Hz. 

The  simulations  were  conducted  using  the  baseband  simulation  model  in  [8], 
operated  at  the  chip  rate.  Thus,  the  sampling  rate  was  set  at  f,  =  Rc,  and  the  decimation 
and  sign  reversal  operations  were  not  required. 

The  phase  and  amplitude  tracking  capabilities  of  the  interference  estimator  are 
shown  in  Figs.  8(a)  and  (b)  for  the  simulation  conditions  in  Table  1.  The  d  profile  for  this 
test  is  in  Fig.  7  and  the  bandwidth  of  the  interference  was  .01  Hz.  Referring  to  Fig.  8(a), 
for  large  values  of  d  the  phase  noise  is  evident  until  at  about  400  iterations,  at  which 
point  the  phase  noise  has  decreased  substantially;  this  indicates  that  the  bandwidth  of 
the  Kalman  filter  has  decreased.  Observe  also  that  initially,  the  algorithm  exhibits  some 
slippage,  taking  from  10  to  15  radians  before  it  settles  down  to  a  tracking  error  of  2t 
radians  after  approximately  50  iterations.  The  initial  slippage  is  due  to  d  being  initially 
too  large.  As  d  becomes  too  small,  the  algorithm  eventually  loses  track  because  the  rate  of 
change  of  the  phase  is  too  fast  for  it  (i.e.,  the  bandwidth  of  the  Kalman  filter  is  too  small). 
Figure  8(b)  illustrates  the  estimate  of  the  envelope  of  the  interference.  As  d  decreases. 


Figure  7:  Profile  of  the  frequency  deviation  constant  d  used  in  the  phase  and  amplitude 
tracking  tests. 
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I{k)  appears  to  improve,  until  at  approximately  1200  iterations,  the  estimate  degrades 
substantially.  The  results  in  Fig.  8  suggest  that  an  optimum  value  for  d  exists  to  achieve 
accurate  tracking  of  the  interferer. 

The  bit  error  rate  as  a  function  of  several  values  of  d  is  illustrated  in  Fig.  9  for  an 
interference  bandwidth  Bi  —  0.005  Hz  and  low  pass  filter  bandwidths  Bipp  =  .10  and  .20 
Hz.  As  an  example,  the  bit  error  rate  decreases  for  increasing  Blpf,  and  the  minimum 
bit  error  rate  shifts  to  the  right,  implying  that  both  the  interference  envelope  and  phase 
estimates  are  more  accurate.  One  can  keep  increasing  Blpf\  however,  at  some  point,  the 
bandwidth  of  the  lowpass  filter  will  be  sufficient  to  provide  a  minimum  bit  error  rate. 
Although  not  shown  in  Fig.  (9),  a  bit  error  curve  was  generated  for  the  case  where  Blpf 
was  set  at  .15  Hz,  but  it  differed  only  slightly  from  the  .20  Hz  curve,  suggesting  that 
the  optimum  bandwidth  for  the  lowpass  filter  was  being  approached.  More  simulations 
would  have  to  be  conducted  to  determine  if  bandwidths  beyond  .20  Hz  would  stairt  to 
deteriorate  the  results. 

Figure  10  illustrates  the  effect  on  bit  error  rate  of  an  interferer  with  larger  band- 
width,  i.e.,  =  .025  Hz,  for  the  same  low  pass  filter  bandwidths  used  in  Fig.  9.  The 

degradation  in  performance  is  quite  evident.  Comparing  Figs.  9  and  10,  for  the  same 
Bipf,  the  required  value  for  d  increases  as  the  interference  bandwidth  increases.  Here 
too,  more  simtilations  would  have  to  be  conducted  to  determine  if  bandwidths  beyond  .20 
Hz  would  continue  to  improve  or  deteriorate  the  results. 

Finally,  Fig.  11  exhibits  the  performance  for  B^pf  =  .20  Hz  and  several  inter¬ 
ference  bandwidths  ranging  from  0.005  to  .025  Hz.  As  increases,  the  bit  error  rate 
increases,  and  the  bandwidth  required  by  the  Kalman  filter  also  increases. 

The  effect  on  the  bit  error  rate  by  modifying  the  carrier-to-noise  ratio  assumed  in 
the  interference  model  (i.e.,  differenct  values  of  or  equivalently  N/2)  for  Blpf  =  .20  Hz 
and  Bi  =  .01  Hz  is  illustrated  in  Fig.  12.  Observe  how  the  point  at  which  the  minimum  bit 
error  rate  occurs  shifts  to  the  left  for  increasing  gun  (i.e.,  decreasing  Nf2).  This  implies 
that  as  the  gain  increases,  the  steady-state  bandwidth  of  the  Kalman  filter  increases  for 
a  given  d,  resulting  in  more  phase  noise.  To  compensate  for  this,  a  lower  value  of  d  is 
required.  The  important  point  to  be  made  here  is  that  one  can  select  particular  values 
for  a  and  <rl  and  vary  d  to  provide  an  adaptive  capability. 

The  bit  error  rate  performance  as  a  function  of  Eb/No  for  several  interference 
bandwidths  and  optimum  values  of  d,  the  latter  of  which  were  obtained  from  curves  similar 
to  those  in  Fig.  9  is  illustrated  in  Fig.  13.  The  bandwidth  Bipf  of  the  low  pass  filter 
was  fixed  at  0.4  Hz.  Observe  that  for  increasing  values  of  H,-,  the  performance  degrades. 
This  degradation  is  due  to  two  things:  (a)  increased  levels  of  phase  noise  because  of  the 
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(a) 


f 


(b) 


Figure  8:  An  example  of  the  tracking  ability  of  the  interference  estimator  using  the  d 
profile  in  Fig.  7  (B,-  =  0.01  Hz):  (a)  Phase;  (b)  Envelope. 
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tion  of  d  and  Blpf  for  Bi  =  0.5%  of  the  chip  rate. 
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Bit  Emr  Rate 


Bit  Elixir  Rate 


Figure  12:  Bit  error  r 
to- noise  ratio,  {CNR 


Eb/NO(dB) 


Figure  13:  Bit  error  rate  as  a  function  of  Eh/No  for  several  values  of  interference  band 
width,  expressed  as  a  percentage  of  the  chip  rate. 


requirement  for  l2U’ger  values  of  d  to  track  the  more  rapidly  varying  interference  (this 
translates  into  more  residual  interference  at  the  output  of  the  interference  canceller  in 
Fig.  1);  and,  (b)  signal  distortion  at  the  output  of  the  interference  canceller. 

6.0  ADAPTIVE  ARCHITECTURES 

To  this  point  an  excisor  using  a  Kalman  filter  has  been  discussed,  with  several  performance 
curves  applied  to  a  general  form  of  interference,  i.e.,  Gaussian  narrowband  interference. 
These  results  were  obtained  with  a  non-adaptive  architecture.  This  section  discusses 
three  possible  structures  which  would  provide  adaptive  capabilities.  The  first  is  more 
appropriate  for  interferers  with  slowly  varying  envelopes;  the  second  is  the  more  general 
structure;  and  the  third  uses  decision  feedback,  given  an  a  priori  data  sequence  to  stabilize 
the  Kalman  parameters. 

6.1  ARCHITECTURE  1  -  CONSTANT  ENVELOPE  INTERFERER 

The  first  architecture  is  illustrated  in  Fig.  14(a),  which  is  an  extension  to  Fig.  1.  The 
main  difference  is  the  inclusion  of  an  envelope  detector  which  provides  a  reference  signal 
for  the  Kalman  filter.  For  large  interference-to-signal-plus-noise  ratio  conditions,  this  is  a 
reasonable  approach. 

How  the  reference  signal  Iref{k)  is  used  by  the  interference  estimator  is  illus¬ 
trated  in  Fig.  14(b).  The  adaptive  algorithm  is  intended  to  work  as  follows.  Given  that  a 
narrowband  interferer  is  present  at  some  unknown  frequency  relative  to  the  spread  spec- 
tnim  signal  carrier,  the  frequency  deviation  constant  d  of  the  Kailman  filter  is  set  at  a 
sufficiently  large  value  so  that  the  algorithm  does  not  lock  on  to  any  signal;  this  value  is 
determined  from  the  parameters  of  the  Kalman  filter  (i.e.,  a  and  al)  and  its  operating 
curves  [8].  The  frequency  deviation  constant  d  is  then  decreased  at  some  rate  using  a 
profile  similar  to  the  one  in  Fig.  7.  Before  the  Kalman  filter  locks  on  to  the  interferer, 
the  error  between  the  reference  /re/(^)  /(^)  will  be  large.  It  will  remain  so  until  d 

approaches  the  vicinity  of  the  optimum  value  in  a  mean-squared  context. 

If  the  error  is  defined  as 


e{k)  =  -  i{k\  (40) 

then  for  a  constant  envelope  interferer,  and  large  interference-to-signal-plus-noise  ratios, 

(t{k)  -  e(k\k  - 1))’ « €{k)/i„,.  (41) 


pn 

Correlator 


(a) 


(b) 


Figure  14:  An  adaptive  architecture  in  which  a  reference  signal  /re/(^)  is  used  by  the 
Kalman  algorithm:  (a)  System  level  diagram;  (b)  Amplitude  estimator,  error  processor 
and  controller. 
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This  can  be  seen  by  noting  &om  Eqs.  (36)  and  (37),  that  the  interference  amplitude  is 

i  =  Icx3s{${k)  —  d(k\k  —  1))  +  noise  terms.  (42) 

By  retaining  the  first  two  terms  of  an  approximation  to  the  cos  term  and  assuming  that 
hej  «  /,  the  desired  result  in  Eq.  (41)  is  obtained.  The  Error  processor  determines 
the  mean  of  Eq.  (41)  while  the  box  labelled  Control  monitors  the  mean-squared-error, 
continually  looking  for  a  minimum.  If  d  gets  too  small,  then  the  mean-squared-error  starts 
to  increase,  implying  that  d  must  reverse  its  direction.  An  example  of  the  performance  of 
this  excisor  for  the  case  of  an  FM  interferer  is  in  [15]. 

6.2  ARCHITECTURE  2  -  GENERAL  NARROWBAND  INTERFERENCE 

This  architecture  is  shown  in  Fig.  15(a).  The  function  which  is  minimized  is  E{e{ky} 
preceding  the  despreader.  This  approach  is  similar  to  minimizing  the  squared  error  func¬ 
tion  resulting  from  the  use  of  a  linear  predictor  excisor  [2].  For  the  linear  predictor  case, 
the  function  that  is  minimized  is 

£{£’(t)}  =  E{Hi(k)  -  i{k))  +  s(k)  +  n(«:)in,  (43) 

where 

i{k)  =  ai{i{k  -  1)  +  s{k  -  1)  +  n{k  -  1))  +  •  ’  • 

+  am(iik  -  m)  -b  s{k  —  m)  -|-  n{k  -  m)).  (44) 

In  Eq.  (44),  the  Cj  are  the  linear  predictor  coefficients  of  the  excisor  and  m  is 
the  order  of  the  linear  predictor.  Using  the  fact  that  s{k)  -|-  n{k)  is  uncorrelated  with 
i{k)  —  i{k),  Eq.  (43)  can  be  rewritten  as 

E{e^ik)}  =  E{{t{k)  -  f(fc))n  +  E{{sik)  +  n{k)n  (45) 

The  objective  then,  is  to  select  the  set  of  coefficients  a,-,  i  =  1,  >  •  •  ,m  in  Eq.  (44)  which 
minimizes  Eq.  (45)  (i.e.,  the  term  ^{(*(l!)  — *(k))*})  using,  for  example,  the  orthogonality 
principle  [16]. 

For  the  Kalman  filter  excisor,  an  equation  like  Eq.  (43)  occurs,  the  main  differ¬ 
ence  being  that  i{k)  is  correlated  with  s{k)  -f  n{k)  since  the  filter  form  of  the  Kalman 
filter  has  been  used.  The  result  will  be  a  degradation  in  performance.  To  circumvent 
this  problem,  one  can  use  the  predictor  form  of  the  Kalman  filter,  which  would  result  in 
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an  equation  like  £)q.  (45).  The  trade-off  here,  however,  is  the  fact  that  the  interference 
estimate  will  be 

t(Jb)  =  I(k  —  1)  coaiuok  -f  —  2))  (46) 

as  illustrated  in  Fig.  15(b).  If  it  is  assumed  that  the  envelope  varies  slowly  enough  relative 

A  A 

to  the  sampling  rate  so  that  I{k  —  1)  A^)t 


t(k)  =  i(k)  coaitJok  +  e{k\k  -  2)). 


(47) 


The  objective,  using  the  Kalman  filtering  approach,  would  be  to  minimize  the  term 
E{(i(k)  —  *(A;))*}  in  Eq.  (45)  by  adjusting  d  and  Blpf  for  fixed  and  a. 

For  architectures  1  and  2,  the  signal-to-noise  ratio  at  the  output  of  the  despreader 
can  be  shown  to  be  [17] 


SNR  = 


_ ^ _ 


(48) 


where 

.  ^i{k)  =  i{k)  -  i{k),  (49) 

and  where  it  has  been  assumed  that  the  filter /prechctor  form  of  the  Kalman  filter  has  been 
used.  With  the  assumption  that  the  noise  at  the  output  of  the  despreader  is  Gaussian, 
the  resulting  BER  is,  therefore. 


Pb  =  lerfcy/S  NR/2, 

A 

where  cr/c(x)  is  the  complementary  error  function. 


(50) 


6.3  ARCHIXECTURE  3  -  DECISION  FEEDBACK 

The  third  architecture  is  illustrated  in  Fig.  16.  The  box  shown  as  Interference  Estimator 
uses  the  filter/predictor  form  (Fig.  15)  without  the  delays  following  bp{k)  and  z{k).  A 
training  signal  is  provided  initially  to  give  time  for  the  Kalman  filter  to  stabilize  without 
the  signal  being  present,  i.e.,  the  input  z{k  —  1)  during  the  training  period  has  the  known 
signal  removed  from  it.  The  output  before  the  limiter  and  bandpass  filter  is,  therefore, 

z\k-l)  =  2(ifc  -  1)  -  s(ib  -  1) 

=  t(fc  —  1) -f- n(fc  —  1).  (51) 
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Ii(k)  iliOt) 


(b)  A 


Figure  15:  An  adaptive  architecture  in  which  £{6^(Jb)}  is  minimized:  (a)  System  level 
diagram;  (b)  Filter/Predictor  form  of  the  Kalman  filter  and  amplitude  estimator. 
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Diiring  the  training  period,  the  interference  estimate  will  be  more  accurate.  The  output 
of  the  pn  correlator,  delayed  by  one  chip,  has  the  known  signal  also  removed  from  it, 
resulting  in  the  error  signal 


e{k-l)=n{k-l)  +  Ai{k-l),  (52) 

where  Ai{k  —  1)  is  the  residual  interference.  The  function  that  is  minizmized  is  E{e{k  — 
1)^}.  When  the  training  sequence  is  completed,  the  system  shifts  to  the  decision  feedback 
mode,  where  the  output  of  the  bit  detector  produces  a  bit  decision  every  Te  seconds, 
which  is  fedback.  This  feedback  signal  is  respread  and  scaled  to  the  chip  energy  level.  It 
should  be  noted  that  the  bit  decision  improves  as  the  bit  evolves  over  the  L  chips  at  the 
output  of  the  summer  following  the  pn  correlation  process;  the  summer  is  cleared  every 
LTc  seconds. 

Once  the  training  sequence  has  passed,  the  performance  will  be  a  function  of  the 
bit  error  rate.  During  this  mode  of  operation,  the  input  to  the  limiter  and  bandpass  filter 
will  be 

z\k  —  1)  =  {s{k  —  1)  —  s{k  —  1))  +  n{k  —  1)  +  i{k  —  1).  (53) 

As  can  be  seen  from  Eq.  (53),  instances  will  occur  (more  so  near  the  beginning  of  each  bit) 
when  the  the  first  term  will  not  be  zero,  resulting  in  some  degradation  in  performance 
compared  to  the  case  when  complete  signal  cancellation  occurs.  Further  work  is  still 
required  to  assess  this  particular  architecture. 


Figure  16:  A  decision-feedback  excisor. 


7.0  CONCLUDEMG  REMARKS 

A  novel  Kalman  filtering  approach  has  been  presented  for  the  suppression  of  interference 
in  direct  sequence  spread  spectrum  communications  systems.  This  approach  differs  from 
other  time-modelling  techniques  based  on  Autoregressive  methods. 

Several  bit  error  rate  curves  were  obt^ned  through  simulation  for  the  case  of 
narrowband  Gaussian  noise  interference.  For  increasing  interference  bandwidths,  it  was 
shown  that  performance  degraded.  It  was  also  shown  that  performance  was  a  function  of 
the  frequency  deviation  constant  d  of  the  Kalman  filter,  for  fixed  filter  gain  and  “message 
bandwidth”  a  of  the  state  space  model.  The  bit  error  rate  was  more  sensitive  to  the 
frequency  deviation  constant  for  quite  narrow  interference  bandwidths  (in  the  results  the 
smallest  bandwidth  used  was  0.005  Hz)  than  for  larger  bandwidths.  For  this  latter  case, 
there  was  a  broad  range  of  values  for  d  which  would  suffice.  It  was  also  shown  that 
performance  was  dependent  on  the  bandwidth  of  the  lowpass  filter  used  in  filtering  the 
envelope  of  the  interference. 

The  bit  error  rate  results  presented  were  not  based  on  2Jiy  particular  adaptive  ar¬ 
chitecture.  The  last  section  of  this  report  presented  three  possible  architectures,  each  one 
minimizing  a  different  squared  error  function.  The  first  architecture  was  more  amenable 
to  interferers  with  a  constant  or  slowly  varying  envelope.  The  second  was  concerned 
with  suppressing  the  more  general  interferer  and  involved  minimizing  the  squared  error 
function  also  used  in  the  AR  techniques.  Finally,  the  third  architecture  incorporated  de¬ 
cision  feedback  coupled  with  an  initial  training  sequence.  The  main  objective  associated 
with  this  last  approach  was  to  reduce  degrading  effects  due  to  the  signal  being  present  in 
estimating  the  interference. 

Two  areas  of  work  remjun  to  be  completed.  These  include  the  development  of 
theoretical  bit  error  rate  curves  for  the  three  architectures,  and  the  assessment  of  the 
decision  feedback  architecture  through  computer  simulation. 
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